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We present a density-functional theory study addressing the on-surface adsorption of oxygen 
at the Pd(llA'') {N = 3,5,7) vicinal surfaces, which exhibit (111) steps and (100) terraces of 
increasing width. We find the binding to be predominantly governed by the local coordination at 
the adsorption site. This leads to very similar bonding properties at the threefold step sites of all 
three vicinal surfaces, while the binding at the central fourfold hollow site in the four atomic row 
terrace of Pd(117) is aheady very little disturbed by the presence of the neighboring steps. 

PACS numbers: 68.43.Bc,71.15.Mb,81.65.Mq 
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I. INTRODUCTION 

Important materials processes and functions like oxide 
formation or oxidation catalysis have motivated a wealth 
of Surface Science studies on low Miller index surfaces of 
late transition metals (TMs) and their interaction with 
oxygeni In recent years, experimental work has been 
increasingly complemented by independent information 
from first-principles theoretical studies, which appropri- 
ately account for the band structure of the TMs through 
supercell setups with periodic boundary conditions liS. En 
route towards an atomic-scale understanding of real TM 
surfaces, investigating the role of defects like steps, kinks 
or vacancies with similar rigor and detail as established 
for the low-index surfaces becomes a key endeavor. With 
respect to catalysis and oxide formation, such defects are 
frequently discussed as playing a prominent role, e.g. in 
form of active sites facilitating molecular dissociation or 
as nucleation centersii*^ 

Steps and kinks are suitably studied at high-index vici- 
nal surfaces, which contain regular arrays of steps of spe- 
cific orientation, separated by terraces of specific width. 
The advantage of these surfaces is that they (ideally) 
only contain one defined type of step in always the same 
local environment, and are amenable to periodic bound- 
ary supercell calculations. The disadvantage is that the 
extracted step properties can be masked by step-step 
interactions at decreasing terrace widths. This holds 
in particular, since one would precisely try to perform 
first-principles calculations addressing a certain step type 
by employing vicinal surfaces with a minimum terrace 
width, in order to reduce the computational burden. In 
this respect, we find it useful to conduct a systematic 
study of a family of vicinal surfaces that always exhibit 
the same step type, but terraces of increasing width. Ex- 
amining the geometric and electronic structure of the 
clean and oxygen-covered surfaces, we can then extract 
by how much the local properties at a given step are af- 
fected by the presence of the neighboring steps. 

Specifically, we use density-functional theory (DFT) to 
study the series of vicinal Pd(llA'^) (7V=3,5,7) surfaces, 
additionally comparing them to the low-index Pd(lOO) 
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FIG. 1: (Color online) Top and side views of the atomic ge- 
ometry of the family of Pd(llA'^) surfaces {N=3 [top panel], 
5 [middle panel], 7 [bottom panel]) surfaces. Light grey (yel- 
low) spheres represent Pd atoms, while the Pd step atoms are 
shown in dark grey (green) , to better illustrate the step struc- 
ture in the top view. Additionally shown in the top view are 
the (1x1) and (1x2) surface unit-cells (solid and dashed lines, 
respectively), as well as the high-symmetry sites for oxygen 
adsorption considered in this study (see text). The side view 
includes also the layer numbering and the relaxed spacings dij 
between layers i and j (and in parenthesis the percent change 
Adij with respect to the bulk interlayer spacing dt). 



and Pd(lll) surfaces. Figure explains the geometric 
structure of the vicinals, which exhibit (111) steps and 
(100) terraces with 2, 3 and 4 atomic rows for Pd(113), 
Pd(115) and Pd(117), respectively. Besides the surface 
energetics and geometric structure with and without oxy- 
gen adsorbates, we also study the local density of states 
(LDOS) and initial-state surface core-level shifts (SCLS) 
as measures of the surface electronic structure. Apart 
from insight into details of the oxygen-metal chemical 
bond, our major conclusion is that the oxygen adsorption 
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properties at the (111) step of Pd(117) (with 4 atomic 
row terraces) are aheady very httlc disturbed by the pres- 
ence of the neighboring steps of the vicinal. 

II. THEORY 

The DFT calculations were performed within the 
highly accurate full-potential augmented plane wave -I- 
local orbitals (LAPW/APW+lo) scheme^, using the gen- 
eralized gradient approximation (GGA-PBE)® for the 
exchange-correlation (xc) functional. All surfaces were 
modeled by supercells, employing symmetric slabs (with 
O adsorption on both sides of the slab) and a minimum 
vacuum thickness of 12 A to ensure the decoupling of 
consecutive slabs. The slab thickness was 7 layers for 
Pd(lll) and Pd(lOO), 9 layers for Pd(113), 17 layers for 
Pd(115) and 19 layers for Pd(119), cf. Fig. □ For the 
oxygen adsorption studies the thickness of the Pd(113), 
Pd(115) and Pd(117) slabs was increased to 13, 17, and 
23, respectively. All structures were fully relaxed, keep- 
ing only the atomic positions in the central two slab layers 
at their fixed bulk positions. 

The LAPW/APW-l-lo basis set parameters are as fol- 
lows: i?MT — 2.1bohr, i?MT — l.lbohr, wave function 
expansion inside the muffin tin spheres up to ^j^^x — 12, 
and potential expansion up to l]^^^ = 6. The energy cut- 
off for the plane wave representation in the interstitial 
region between the mufhn tin spheres was E'^^^ = 20 Ry 
for the wave functions and E^^^ = 196 Ry for the po- 
tential. The core states are treated fully relativistic, 
while the semi-core and valence states are treated within 
a scalar relativistic approximation, i.e. spin-orbit cou- 
pling is included (neglected) for the core (semi-core and 
valence) states. Monkhorst-Pack (MP) grids were used 
for the Brillouin zone integrations. Specifically, we used 
(6 X 10 X 1), (4 X 10 X 1) and (3 x 10 x 1) grids for the cal- 
culations of (1 X 1) surface unit-cells of Pd(113), Pd(115), 
and Pd(117), respectively. The selected calculations on 
Pd(lOO) and Pd(lll) used for the comparison with low- 
index surfaces were done with a (12x12x1) grid. For the 
larger surface cells, the grids were reduced accordingly to 
keep the same sampling of the reciprocal space. 

The central energetic quantities evaluated in our study 
are the surface energy 7 for the clean surfaces, and the 
binding energy Et for the oxygen adsorption. The prior 
is defined as 

7 = ^ (£^slab - NsEhnlk) , (1) 

where -Esiab and -Ebuik are the total energy of the slab and 
of a bulk atom, respectively. is the number of atoms 
in the slab in the employed supercell, A the correspond- 
ing surface area, and the factor of ^ is used to account 
for the two equivalent surfaces of the slabs. In order to 
minimize the error due to the finite k-point sampling, 
^'buik is obtained for each surface from a calculation in 
exactly the same supercell as used for -Esiab, and filling 



the vacuum region with additional Pd layers to generate 
a bulk system. The oxygen binding energy is defined as 

Eb - [eo/p^-Ep^-^Eo,^ , (2) 

where A'^o is the number of O atoms present in the em- 
ployed supercell, and Eq/p^, Ep^ and E02 are the total 
energies of the slab containing the oxygen adsorbate, the 
corresponding clean Pd slab, and the total energy of an 
isolated molecule, respectively. Again, the factor i ac- 
counts for the fact that we adsorb oxygen at both sides 
of the symmetric slab. 

With this definition, a positive Ei, value indicates that 
dissociative adsorption of O2 is exothermic at T = OK. 
To obtain the total energy of the isolated O2 molecule, 
we exploit the relation Eq^ = 2E^°"' - D, where 
is the total energy of an isolated oxygen atom, and D the 
theoretical O2 binding energjiS. The isolated O atom is 
then calculated spin-polarized, inside a rectangular cell 
of side lengths 12 x 13 x 14bohr, F-point sampling of 
the Brillouin zone and without spherically averaging the 
electron density in the open valence shell. For D we 
employ the previously published ultra-converged value of 
6.21 eV^. Compared to the experimental binding energy 
of 5.12 eViSi, this value shows the well-known, substantial 
GGA-induced error. We note, however, that this will 
neither affect our discussion, nor our conclusions, since 
we focus on differences between oxygen binding energies 
at different surfaces and different sites, where this error 
cancels out. 

Limitations in the numerical accuracy of our results 
come predominantly from the finite LAPW/APW-l-lo ba- 
sis set and the supercell geometry. We performed ex- 
tensive test calculations to estimate the importance of 
the various numerical approximations. Our tests concen- 
trated on the convergence of the crucial energetic quan- 
tities of our study, using the surface energies of all five 
surfaces, as well as the oxygen binding energy in the Thu 
and Sh2 site (see definition below) in (1 x 1) unit-cells 
on the vicinals as representatives. Naming here only the 
factors that are most critical for the convergence, we in- 
creased the employed cutoff for the plane-wave basis set 
systematically from 20 to 30 Ry, doubled the k-mesh den- 
sity, and tested slab thicknesses with up to twice as many 
layers. From all these tests we estimate that our reported 
surface energies are converged to within 2 me V/A2, and 
the O binding energies to within 150mcV/O atom. How- 
ever, much more important for our study comparing the 
properties at different surfaces and sites are relative bind- 
ing energy and surface energy differences. For these dif- 
ferences, a large part of the uncertainty, e.g. due to 
the description of the isolated O atom, cancels out, and 
our tests show that the numerical uncertainty for Aii'f, 
of O atoms at different sites or at different surfaces is 
±50meV/O atom, and the relative surface energy differ- 
ences are converged to within 0.5 me V/A2. As will be- 
come apparent below, this does not affect our discussion 
and conclusions. 
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III. CLEAN VICINAL SURFACES 

A. Geometric structure 

Using the Murnaghan equation of state, we com- 
pute the DFT-GGA-PBE lattice constant for Pd as 
a ~ 3.947 A (neglecting zero-point vibrations) and the 
bulk modulus as i? = 157 GPa. These results are in 
excellent agreement with our previous study using the 
earlier WIEN97 LAPW code, which gave a = 3.944 A 
and B = 163 GPa, respectiveljiii. The slight overes- 
timation of the lattice constant (about +2% compared 
to the experimental value of Ooxp = 3.89 Ai^) and the 
(corresponding) underestimation of the bulk modulus 
{Bcxp = 181 GPaiSi) are in line with analogous studies 
for other late Ad TMsiiSiiii^ The obtained surface re- 
laxation pattern for the low-index Pd(lOO) and Pd(lll) 
surfaces is also as expected in view of the trend under- 
standing over the late 4d TMs e.g. discussed by Meth- 
fessel, Hennig and Schcfllcr — . For both surfaces we ob- 
tain only insignificant (< 1 %) variations in the interlayer 
spacings of the topmost layers with respect to the bulk 
value. 

Turning to the surface relaxation for the three vicinal 
surfaces, our results are summarized aside with the ge- 
ometry side views in Fig. ^ One can clearly discern that 
the largest deviations from the bulk spacing are concen- 
trated to those layers, in which the atoms have a lower 
coordination than in the bulk. When performing our 
slab thickness tests, we also fully relaxed slabs with up 
to twice as many layers as those shown in Fig. ^ How- 
ever, no changes in the interlayer spacings were observed 
compared to the values listed in Fig. ^ and the deeper 
layers showed insignificant relaxations. Focusing there- 
fore on the larger relaxations in the topmost layers, we 
notice that the signs of the percent changes Ad^ of the 
interlayer spacings compared to the bulk spacing exhibit 
a certain pattern. Starting with the topmost layer dis- 
tance, we namely have — I — for Pd(113), 1 — for 

Pd(115) and h - for Pd(117), where "-" corre- 
sponds to a contraction and "+" to an expansion of the 
interlayer spacing. 

There are two theoretical models that explain these 
sequences as a consequence of Smoluchowski's ideas of 
charge smoothing (or equivalently, in a chemical lan- 
guage, through the bond-order bond-length correlation 
concept )iSiiLi2>: The atom rows modeii^ predicts that for 
a surface with n^ow atomic rows on the terraces, the first 
(rirow — 1) topmost interlayer spacings exhibit a contrac- 
tion, while the riiowth one expands. The nearest-neighbor 
(NN) coordination model^^ focuses instead on those sur- 
face layers which have a NN coordination smaller than 
the bulk. For a surface where the tinn topmost layers 
are undercoordinated, the topmost (tinn — 1) interlayer 
spacings are then expected to contract, while the riNNth 
one expands. As apparent from Table ^ the multilayer 
relaxation sequences obtained for the present family of 
Pd(llA^) vicinals conforms with both models, which is 



TABLE I: Comparison of the multilayer relaxation sequence 
of the Pd(117V), (A'' = 3,5,7) vicinals to the key quantities 
of the atom-rows model— and NN coordination model^°. A 
"— " in the relaxation sequence stands for a contraction of the 
corresponding interlayer spacing, and a for an expansion, 
nrow is the number of atomic rows on the terraces and unn 
the number of undercoordinated layers at the surface. For 
comparison, the coordination is 12 in the fee bulk, and 8 for 
a (100) terrace atom. 



Vicinal 


TT-row 


NN Coordination 


riNN 


Relaxation 


Surface 




For Topmost Layers 




Sequence 


Pd(113) 


2 


7, 10, 12, 12, ... 


2 


- + - 


Pd(115) 


3 


7, 8, 10, 12, 12, ... 


3 





Pd(117) 


4 


7, 8, 8, 10, 12, 12, ... 


4 


h - 



TABLE IL Computed surface energies of the low-index 
Pd(lOO) and Pd(lll) surfaces, as well as of the three vicinal 
(Pd(lliV), {N = 3,5,7)) surfaces. Bulk-truncated refers to 
a bulk-truncated surface geometry, while relaxed corresponds 
to a fully relaxed surface geometry. We present all values in 
meV/A^ and in eV/surface atom. 



meV/A^ Pd(lll) 


Pd(113) Pd(115) Pd(117) 


Pd(lOO) 


Bulk-truncated 87.9 
Relaxed 87.9 


100.3 99.4 99.1 
98.3 97.3 97.4 


96.4 
96.3 


eV/atom Pd(lll) 


Pd(113) Pd(115) Pd(117) 


Pd(lOO) 


Bulk-truncated 0.59 
Relaxed 0.59 


1.30 2.01 2.76 
1.27 1.97 2.71 


0.75 
0.75 



no surprise, since for the case of (111) steps the atom- 
rows model is actually a special case of the more general 
NN coordination concept. 



B. Energetics and electronic structure 

Table Hll summarizes the computed surface energies of 
the three vicinal surfaces and compares them to the sur- 
face energies of the two low-index surfaces Pd(lOO) and 
Pd(lll). Due to the negligible layer relaxation of the 
latter two surfaces, a full relaxation of the surface ge- 
ometry has only a negligible effect on their surface en- 
ergy. As expected, this is different for the more open 
vicinal surfaces, where the atomic relaxation leads to 
some lowering. The values obtained for Pd(lll) and 
Pd(lOO), 0.59 eV/atom and 0.75 eV/atom respectively, 
are in good agreement (±5 %) with recent other DFT- 
GGA studiesi^iSLS^. However, they are significantly dif- 
ferent to two earlier trend studies, namely the linear- 
mufhn-tin-orbital (LMTO) GGA work by Vitos et al. 
(Pd(lll) 0.82 eV/atom; Pd(lOO) 1.15 eV/atom)2i, and 
the LMTO work by Methfessel, Hennig and Scheffler 
employing the local-density approximation (LDA)-- to 
the xc functional (Pd(lll) 0.68 eV/atom; Pd(lOO) 0.89 
eV/atom)i^. To assess the effect of the different approx- 
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FIG. 2: Local density of states (LDOS) in the different layers 
in Pd(117). Compared to the bulk-like LDOS in the central 
layer of the slab, the topmost four surface layers (correspond- 
ing to undercoordinated step and terrace atoms) exhibit a 
significant narrowing of the valence 4(i band. The ensuing 
shift of the center of gravity of the valence 4d band induces a 
non-negligible shift Af^ij^fj in the 3d core-level positions with 
respect to the value in a bulk Pd atom. Note the similarity 
of the Afj^^fi exhibited by the two terrace atoms in Pd(117), 
2nd and 3rd row, to the one of the topmost layer atoms in an 
infinite Pd(lOO) terrace, Af.^t'i'fi = -0.50 eV. 



imate xc functional in the latter study, we recomputed 
the surface energies with the LDA and using the opti- 
mized LDA lattice constant (3.84 A). The obtained val- 
ues of 0.77eV/atom for Pd(lll) and 1.02eV/atom for 
Pd(lOO) are again in good agreement with recent other 
LDA studies^'*'^- . However, although the values are sig- 
nificantly different to the GGA numbers, they cannot 
resolve the discrepancy with respect to the early LDA 
work by Methfessel and coworkers. Instead, it seems to 
be the non-relativistic treatment of the valence and semi- 
core electrons in the latter work that is behind the dif- 
ference. When deliberately switching off the scalar rela- 
tivistic treatment employed in our work, we obtain values 
of 0.69eV/atom (Pd(lll)) and 0.88eV/atom (Pd(lOO)) 
that are now very close to those obtained by Methfessel 
and coworkers. While we are therefore able to resolve 
the difference with respect to this work, we are unable 
to account for the discrepancy with respect to the TB- 
LMTO work by Vitos et at. As already discussed by 
Da Silva, Stampfl and Schefllepi^, the reason might be 
that the atomic sphere approximation (ASA) employed 
by Vitos et al. does not allow for sufficient flexibility of 
electronic relaxation at the surface. 



TABLE IIL InitiaLstate surface core-level shifts Afnitifi for 
the topmost layers with reduced NN coordination, cf. Table 
Q at the low-index and vicinal surfaces. All values are given 
in eV. 



A SCLS 
^initial 


Pd(lll) 


Pd(113) 


Pd(115) 


Pd(117) 


Pd(lOO) 


1st 


-0.39 


-0.56 


-0.50 


-0.54 


-0.50 


2nd 




-0.19 


-0.43 


-0.45 




3rd 






-0.23 


-0.48 




4th 








-0.24 





As for the vicinal surfaces, the sequence of Pd(113), 
Pd(115) and Pd(117) surface energies gives already a first 
hint as to the influence of neighboring steps. In this re- 
spect, the obtained rather similar values in me V/A2, cL 
Table suggest that step-step interactions are rather 
small despite the small terrace widths of these surfaces. 
This finding is in full agreement with a preceding tight- 
binding study by Raouafi et ai, which also determined 
neghgible step-step interactions at p(lOO) x (111) vicinals 
of Rh, Pd and CiiS&. The fast decay of the local perturba- 
tion introduced by the (111) step is also nicely discerned 
in the surface electronic structure. Figure[21shows the lo- 
cal density of states (LDOS) calculated inside the muffin 
tin spheres of the atoms in the various surface layers at 
Pd(117). Compared to the bulk-like LDOS of the central 
layer in the slab, only the topmost four layers correspond- 
ing to step or terrace atoms with reduced NN coordina- 
tion show a significant variation, predominantly in form 
of a narrowing of the valence Ad band. In order to main- 
tain local charge neutrality, the center of gravity of the 
more than half full d band has correspondingly moved 
up compared to the situation in a bulk Pd atom, which 
went hand in hand with a repulsive contribution to the 
Kohn-Sham potential^. This potential change is felt by 
the core electrons, as well, and induces a shift in the posi- 
tion of the core-levels. We define this initial-state surface 
core-level shift (SCLS) with respect to a bulk atom as 



A 



SCLS 
initial 



rfacc 



^bulkl 



(3) 



where g^"'"^'^'^® and e]^^^^ are the Kohn-Sham eigenvalues of 
the core state c in a surface and a bulk atom. With this 
definition, the valence d band narrowing at the surface 
leads to a negative shift, and we focus in this work in 
particular on the SCLS exhibited by the 3d Pd levels. 

^hStiai contribution to the measurable SCLS in 

x-ray photoemission spectroscopy (XPS), which in ad- 
dition comprises also the screening contribution of the 
valence electrons in response to the created core hole2i. 
However, it is not so much this connection to an experi- 
mentally accessible quantity that makes us interested in 
the initial-state SCLS, but rather that the Afj^^^'j'fj are 
a very sensitive probe of changes in the local electronic 
structure of an atom in different environments. In this re- 
spect we note that the A^^^^'j'^j exhibited by the eightfold 
coordinated 2nd and 3rd layer terrace atoms at Pd(117), 
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TABLE IV: Binding energies of oxygen atoms at the dif- 
ferent highly-coordinated terrace and step sites of the three 
Pd(llAi') {N=3,5,7) surfaces, see Fig. 0for an illustration of 
all sites and the employed nomenclature. Listed are results 
for the adsorption of one oxygen atom per (1 x 1) and (1 x 2) 
surface unit-cell. At the Pd(113) surface, the Thu and Thl 
sites coincide, which is why only one value is given. All values 
are in eV. 



(1x1) 


Shi 


Sh2 


Thu 


Thl 


The 


Pd(113) 


0.82 


0.87 


0.60 






Pd(115) 


0.80 


0.89 


0.84 


0.52 




Pd(117) 


0.83 


0.94 


0.86 


0.57 


0.80 


(1x2) 


Shi 


Sh2 


Thu 


Thl 


The 


Pd(113) 


1.08 


1.17 


0.96 






Pd(115) 


1.04 


1.17 


1.27 


0.92 




Pd(117) 


1.05 


1.22 


1.27 


0.98 


1.25 



cf. Fig. 121 are already quite similar to the one we com- 
pute for a topmost layer atom in an infinite Pd(lOO) 
terrace, — 0.50eV/atom. In contrast, all deeper layers 
with 12fold bulk-like coordination exhibit only insignifi- 
cant shifts. Table UTTl shows that the nearest-neighbor co- 
ordination is in fact the overriding parameter determin- 
ing the initial-state SCLSs at all surfaces studied here. 
Comparing to the NN coordination listed in Tabled one 
even finds a roughly linear relation of the A?J?^'t'^j with 
the number of NN Pd atoms, which indicates that per- 
turbations further away than the immediate NN shell are 
rapidly screened away. 

IV. OXYGEN ADSORPTION AT PD(lliV) 
VICINAL SURFACES 

A. Binding energy 

Interested in evaluating the effect of the nearby steps 
in the series of Pd(lliV) surfaces, we consider the ad- 
sorption of oxygen as a probe of the local chemical prop- 
erties. For this we focus on the binding of atomic O 
to the highly-coordinated high-symmetry sites offered by 
the vicinal surfaces, namely the fourfold coordinated hol- 
low sites on the terraces and the threefold coordinated 
hollow sites at the steps, cf. Fig. ^ This yields two 
step sites at each vicinal surface, the Shi site where the 
oxygen atom is coordinated to one step atom and the 
Sh2 site where the oxygen atom is coordinated to two 
step atoms. At the Pd(117) surface, there are three ter- 
race hollow sites, one at the lower edge of the step (Thl), 
one at the upper edge of the step (Thu) and one in the 
center of the terrace (The). Due to the reduced terrace 
width, there are only Thl and Thu sites at Pd(115), and 
at Pd(113) the lower and upper edge terrace hollow sites 
coincide. To analyze the effect of coverage we calculated 
the binding energy of one oxygen atom at any of these 



sites both in (1 x 1) and (1x2) surface unit cells. In 
the prior, all sites of one type are then covered, e.g. all 
Shi sites along a step, while in the latter occupied sites 
alternate with empty ones, e.g. one Shi site covered, the 
neighboring one empty etc. 

The obtained binding energies are compiled in Table 
IIVI and can be understood in a rather simple and lo- 
cal picture. Central for the stability at a given site are 
the number of directly coordinated surface Pd atoms and 
their respective Pd coordination (nfold coordinated ter- 
race or step atoms). Sites which are identical with respect 
to these criteria exhibit binding energies that are almost 
degenerate within our numerical accuracy (±50meV/O 
atom). Prominent examples are the similar binding en- 
ergies at the Shi site at all three surfaces (binding to 
one Tfold coordinated step and two lOfold coordinated 
terrace atoms) and at the Sh2 site at all three surfaces 
(binding to two Tfold coordinated step and one lOfold 
coordinated terrace atom). Due to the narrow terrace 
widths this equivalence is only established for the Thu 
(and separately for the Thl) sites at the Pd(115) and 
Pd(117) surfaces. There, Thl corresponds to binding to 
two 8fold coordinated and two lOfold coordinated ter- 
race atoms, while Thu corresponds to binding to two 
Tfold coordinated step atoms and two 8fold coordinated 
terrace atoms. The binding site The in the center of 
a terrace (coordination to four 8fold terrace atoms) de- 
velops only at the Pd(llT) surface, but is then already 
comparable to the fourfold hollow site at Pd(lOO), corre- 
sponding to an infinite terrace width. We arrive at this 
assessment from calculating (3x1) and (4x1) over layers 
at Pd(lOO) with one oxygen atom per surface unit-cell. 
The obtained binding energies are 0.84 eV/0 atom and 
0.85 eV/0 atom, respectively, which shows that the lat- 
eral interactions between the neighboring O atoms along 
the long axis of the unit-cell are already negligible. On 
the other hand, these binding energies are very close to 
the one computed for the The site at Pd(llT), 0.80 eV/0 
atom, revealing an only marginal influence of the neigh- 
boring steps in the latter case. Furthermore, the pre- 
dominant role of the local coordination seems to hold 
independently of the coverage, since our calculations in 
(1x2) surface unit-cells, cf. Table Hvl show exactly the 
same pattern. In this case, we compute an O binding en- 
ergy of 1.29 eV/0 atom in a (3 x 2) overlayer on Pd(lOO), 
which is again very similar to the value of 1.25 eV/0 atom 
at the The site. 

Also the stability ordering among the different site 
types largely follows a local coordination picture, when 
recognizing that the bonding is stronger to lower- 
coordinated Pd atoms than to higher-coordinated ones. 
Correspondingly, the Sh2 site is more stable than the 
Shi (bonding to two step atoms vs. bonding to one step 
atom), while the ordering among the terrace sites goes 
with decreasing stability as Thu, The and Thl. Interest- 
ingly, the most stable terrace site (Thu) is slightly more 
stable than the most stable step site (Sh2) at the lower 
(1x2) coverage studied, while this order is reversed at the 



higher (1 x 1) coverage. This might again be explained in 
a local coordination picture, when realizing that in the 
(1x1) arrangement, neighboring O atoms in rows of any 
of the terrace sites share two Pd atoms, while in rows of 
either of the two step sites they share only one Pd atom. 
In terms of a bond-order conservation picture, it is then 
intelligible that equivalent terrace sites are ~ 0.4eV/O 
atom less stable in the denser (1x1) over layer compared 
to the (1x2) one, while for the step sites this reduc- 
tion is only ~ 0.25 eV/0 atom. With the Thu site only 
slightly more favorable than the Sh2 site in the (1 x 2) 
arrangement, the smaller reduction for the Sh2 site when 
going to the dense (1x1) arrangement leads then to 
slight preference for the Sh2 site. We speculate that this 
change in the stability ordering of terrace and step sites 
with coverage may lead to interesting ordering phenom- 
ena. Addressing this topic requires a proper evaluation 
of the partition function, on which we focus in ongoing 
work. 



B. Geometric and electronic structure 

The calculated relaxed surface geometries reflect the 
local coordination picture developed in the previous sec- 
tion. Focusing first on the 0-Pd bond lengths, the data 
for the (1 X 1) O overlayers compiled in Table IVl clear Iv 
reveals two trends (equivalent results are obtained for 
the (1 X 2) O overlayers): Longer 0-Pd bond lengths 
are found at higher coordinated adsorption sites (fourfold 
terrace sites vs. threefold step sites), and shorter bond 
lengths result in bonds to lower coordinated Pd atoms 
(7fold coordinated step atoms vs. 8 or lOfold coordi- 
nated terrace atoms). Equivalent sites at the different 
surfaces that possess the same coordination and hence 
similar binding energies as discussed above, also exhibit 
similar bond lengths. The local geometry at the Shi 
and Sh2 sites is thus the same at all three vicinials, and 
so is the one at the Thu and Thl sites at Pd(115) and 
Pd(117). Since the bonds to lower coordinated Pd atoms 
are shorter, the threefold step sites Shi and Sh2 are not 
exactly threefold, but a slightly shorter bond length re- 
sults to the step atoms. The same occurs for the terrace 
sites, where the inequivalence of the Thl site at the bot- 
tom of the step is most pronounced: Here, the O atom is 
significantly shifted away from ideal fourfold site, away 
from the step and towards the terrace center. Just as 
with respect to the binding energies, the The site is al- 
ready very much comparable in its local geometry to the 
fourfold hollow site at a similarly covered Pd(lOO) sur- 
face: For (3 X 1) and (4 x 1) overlayers with one O atom 
per surface unit-cell we compute bondlengths of 2.14 A, 
which are almost identical to the 2.13A/2.16A found at 
Pd(117). Interestingly, this bond length similarity to the 
low-index surfaces holds even for the threefold step sites 
Shi and Sh2, which exhibit bond lengths that are very 
close to the « 2 A computed for the threefold hollow sites 
in various overlayers on Pd(lll)ii. 
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FIG. 3: (Color online) Geometry side views for adsorption 
in the Sh2 (left panels) and Thu (right panels) sites in (1 x 1) 
overlayers at the Pd(113) [top panels], Pd(115) [middle pan- 
els], and Pd(117) [bottom panels]. The insets show the per- 
cent change of the layer spacings with respect to the bulk 
spacing Adij. For comparison, the value in parenthesis gives 
the corresponding Ad^j of the clean surface. The arrows in- 
dicate the direction of the substrate atom displacement (in- 
cluding the lateral relaxation) with respect to the position in 
the clean surface. Light grey (yellow) spheres represent Pd 
atoms, dark grey (green) spheres Pd step atoms, and small 
dark (red) spheres O atoms. 



The chemisorption of oxygen also induces substantial 
changes in the relaxation pattern of the underlying Pd 
substrate atoms. Similar to the situation at lower cov- 
erages at the low-index Pd(lll)^^ and Pd(lOO) surfaces, 
the main effect of the adsorbed oxygen is to pull the di- 
rectly coordinated Pd atoms out of the surface. However, 
while primarily at Pd(lOO) this leads simply to a signifi- 
cant expansion of the first layer spacing (A(ii2 = +12% 
at (1 X 1) coverage), the relaxation is more complicated 
at the more open vicinal surfaces, which also allow for 
quite some lateral displacements of the Pd atoms. The 
corresponding data is shown in Fig. 13 for O adsorption 
in the most stable step and terrace site, Sh2 and Thu 
respectively, and compiled in Table IVII for all adsorp- 
tion sites. Due to the lateral displacements, the oxygen- 
induced changes in the relaxation pattern also propagate 
much further than just the nearest-neighbor Pd shell. 
Here, the locality is thus not that pronounced, but ap- 
parently this has no significant consequences on the afore 
discussed binding energies. 

The perturbation in the electronic system created by 
the chemisorbed oxygen atom is again very localized, as 
is nicely reflected in the layer-resolved LDOSs. Figure 
0] exemplifies this for oxygen at the Sh2 step and The 
terrace site at Pd(117). In both cases, the local density 
of states is primarily affected in the surface layers cor- 
responding to the Pd atoms directly coordinated to the 
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TABLE V: O-Pd bond lengths (in A) at the different sites in (1 x 1) overlayers. The shorter bond lengths are always to the 
lower coordinated Pd atoms, cf. Fig. 0and Table ITTl At the Pd(113) surface, the Thu and Thl sites coincide, which is why 
only one set of values is given. 





Shi 


Sh2 


Thu Thl 


The 


Pd(113) 
Pd(115) 
Pd(117) 


1.97/2.02/2.02 
1.99/2.01/2.01 
1.98/2.01/2.01 


1.98/1.98/2.02 
1.99/1.99/2.03 
1.98/1.98/2.03 


2.00/2.00/2.47/2.47 
2.15/2.15/2.14/2.14 2.05/2.05/2.30/2.30 
2.15/2.15/2.14/2.14 2.05/2.05/2.30/2.30 


2.13/2.13/2.16/2.16 



TABLE VI: Surface relaxation pattern of the Pd(117V) (A'^ = 3, 5, 7) vicinal surfaces, both for the clean surfaces and with a 
(1 X 1) O overlayer adsorbed. Shown are the interlayer spacings dij (see Fig. Qfor a definition of the layer numbers), as well 
as in parenthesis the percent change Adij with respect to the interlayer spacing di, in the bulk. See Fig. 0for an illustration 
of the various adsorption sites. 



Pd(113) 



Shi 



Sh2 



1.09 A (-9%) 
1.32 A (+11%) 
1.11 A (-7%) 
1.23 A (+3%) 
1.17 A (-2%) 



1.15 A (-4%) 
1.31 A (+10%) 
1.15 A (-3%) 
1.21 A (+2%) 
1.19 A (0%) 



Thu/Thl 



clean 



1.06 A (-11%) 
1.25 A (+5%) 

1.18 A (-1%) 

1.19 A (0%) 
1.19 A (0%) 



(Adas) 

(Ad34) 
(Ad45) 

(Adse) 



1.10 A (-8%) 
1.32 A (+11%) 
1.12 A (-5%) 
1.23 A (+3%) 
1.18A (-1%) 



Pd(115) 


Shi 


Sh2 


Thu 


Thl 


clean 


di2 (Adi2) 

d23 (Ad23) 
d34 (Ad34) 
d45 (Ad45) 

dse (Adse) 
d67 (Ade?) 
drs (Adrs) 


0.74 A (-2%) 
0.62 A (-19%) 
0.95 A (+24%) 
0.70 A (-7%) 
0.71 A (-6%) 
0.83 A (+9%) 
0.72 A (-5%) 


0.79 A (+5%) 
0.58 A (-24%) 
0.97 A (+27%) 
0.71 A (-6%) 
0.72 A (-5%) 
0.81 A (+7%) 
0.74 A (-2%) 


0.60 A (-21%) 
0.71 A (-6%) 
1.00 A (+31%) 
0.55 A (-27%) 
0.85 A (+12%) 
0.77A (+1%) 
0.74 A (-2%) 


0.64 A (-15%) 
0.73 A (-4%) 
0.82 A (+8%) 
0.85 A (+11%) 
0.62 A (-19%) 
0.81 A (+6%) 
0.80 A (+5%) 


0.63 A (-17%) 
0.63 A (-17%) 
0.93 A (+22%) 
0.68 A (-11%) 
0.74 A (-2%) 
0.79 A (+4%) 
0.75 A (-2%) 



Pd(117) 


Shi 


Sh2 


Thu 


The 


Thl 


clean 


di2 (Adi2) 

d23 (Ad23) 
d34 (Ad34) 
d45 (Ad45) 

dse (Adse) 
der (Ads?) 
drs (Adrs) 

dgg (Adgg) 


0.66 A (+20%) 
0.52 A (-6%) 
0.29 A (-47%) 
0.80 A (+46%) 
0.59 A (+8%) 
0.45 A (-18%) 
0.49 A (-11%) 
0.69 A (+25%) 


0.66 A (+20%) 
0.54 A (-1%) 
0.35 A (-37%) 
0.75 A (+36%) 
0.57 A (+4%) 
0.50 A (-9%) 
0.52 A (-6%) 
0.61 A (+12%) 


0.44 A (-21%) 
0.57 A (+3%) 
0.46 A (-17%) 
0.78 A (+42%) 
0.39 A (-30%) 
0.58 A (+5%) 
0.58 A (+5%) 
0.58 A (+5%) 


0.48 A (-13%) 
0.51 A (-7%) 
0.51 A (-7%) 
0.65 A (+18%) 
0.64 A (+17%) 
0.36 A (-34%) 
0.63 A (+15%) 
0.59 A (+6%) 


0.48 A (-12%) 
0.46 A (-16%) 
0.51 A (-8%) 
0.70 A (+28%) 
0.49 A (-11%) 
0.66 A (+20%) 
0.40 A (-28%) 
0.60 A (+10%) 


0.49 A (-12%) 
0.50 A (-10%) 
0.49 A (-12%) 
0.66 A (+20%) 
0.52 A (-6%) 
0.52 A (-6%) 
0.56 A (+1%) 
0.56 A (+1%) 



O adsorbate, namely layers 1 and 4 for the Sh2 site, and 
layers 2 and 3 for the The site, cf. Fig. ^ In these layers, 
we find 0-Pd bonding states just below and antibonding 
states just above the nearly fully occupied Pd Ad band. 
In the shown LDOS for the (1 x 1) overlayer, the bond- 
ing states show some substructure due the presence of a 
noticeable oxygen-oxygen interaction along the rows in 
the short direction of the surface unit-cell. This leads 
to the formation of an adsorbate band structure as pre- 
viously discussed for high-coverage O overlayers at the 
low-index Pd(lll) surfaceii. The surface layers corre- 
sponding to atoms not belonging to the immediate NN 
shell of the oxygen adsorbate show only insignificant vari- 
ations in the LDOS compared to the case of the clean sur- 
face, as does the LDOS in the deeper layers. Correspond- 
ingly, it is also only the directly coordinated surface lay- 
ers that exhibit appreciable adsorbate-induced changes 



in the initial-state SCLSs, cf. Fig. 01 The widening of 
the d-band leads then to more positive shifts compared 
to the clean surface SCLSs^l. The stronger the formed 
0-Pd bond, the larger the disturbance of the local po- 
tential and the larger the ensuing shift. In this respect, 
we find 0-induced shifts, i.e. a change of A^^^^'j'^i com- 
pared to the value at the clean surface, of +1.21 eV and 
+0.48 eV at the 7fold and lOfold coordinated Pd atoms 
to which the O atom binds in the Sh2 site, and shifts of 
+0.69 eV and +0.82 eV at the two types of 8fold terrace 
Pd atoms to which the O atom binds in the The site. This 
indicates a significantly different bonding of the O atom 
to the different NN Pd atoms, which correlates with the 
varying 0-Pd bond lengths shown in Table and sug- 
gests that the 7fold coordinated step atoms may offer 
a stronger bonding compared to the higher-coordinated 
terrace atoms. That the overall binding energy at the 
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FIG. 4: (Color online) Oxygen-induced changes in the local density of states (LDOS) and initial-state SCLSs. Shown are 
the data for oxygen adsorption into the Sh2 step (left panels) and The terrace site (right panels) at Pd(117). Dashed lines 
correspond to the LDOS with the O adsorbate (inside the O muffin tin in the top panel and inside the Pd muffin tin in all 
other panels), solid lines correspond to the clean surface. The Af^yfi exhibited by the Pd terrace atoms directly coordinated 
to the O atom in the The site, layers 2 and 3, are very similar to the 0.32 eV computed for the directly-coordinated Pd atoms 
in a (3 X I) O-overlayer at Pd(IOO). The values in parenthesis give the initial-state SCLSs at the clean surface for comparison. 



Sh2 site is nevertheless still comparable to the most sta- 
ble terrace sites, cf. Table Hvl results simply because the 
latter offer a higher (fourfold) coordination to Pd atoms 
compared to the threefold sites at the (111) step. 

The perturbation due to oxygen adsorption at all other 
investigated sites is as localized as in the just discussed 
case of the Sh2 and The sites. In all cases, we only find 
significant 0-induced shifts in the initial-state SCLS at 
the directly coordinated Pd atoms. Just as for the bind- 
ing energies and 0-Pd bond lengths, also these SCLSs 
are very similar for adsorption into equivalent sites at the 
three vicinal surfaces. In a (3 x 1) overlayer at Pd(lOO) 
we compute an 0-induced SCLS shift of -f0.82eV exhib- 
ited by the directly coordinated Pd atoms, which is again 
very similar to the above cited values for O adsorption 
into the The terrace site at a Pd(117) surface. This cor- 
roborates our understanding that the neighboring steps 
have already a negligible influence on the O bonding at 
this site in the center of an only four atomic row wide 
terrace. 



V. CONCLUSIONS 

We presented a DFT-GGA study addressing the on- 
surface adsorption of oxygen at Pd(llAf) (iV=:3,5,7) vic- 
inal surfaces. Apart from providing detailed insight into 
the oxygen-metal chemical bond, our particular interest 
was to investigate the influence on the adsorption prop- 



erties by the neighboring steps that are separated by 2, 
3 and 4 atomic row wide terraces at the three surfaces, 
respectively. The computed data on binding energies, 
local density of states and initial-state surface core-level 
shifts points at a very localized perturbation created by 
the chemisorbed oxygen atoms, primarily concentrating 
on the nearest neighbor Pd atoms. Central for the sta- 
bility at a given site are then the number of directly co- 
ordinated surface Pd atoms (threefold hollow sites at the 
(111) steps vs. fourfold hollow sites at the (100) terraces) 
and their respective Pd coordination (nfold coordinated 
Pd terrace or step atoms). Sites which are equivalent 
with respect to these criteria exhibit binding energies 
that are degenerate to within 70meV/O atom. Exam- 
ples are the step sites at all three vicinal surfaces, while 
the hollow site in the center of the four atomic row (100) 
terrace at Pd(117) is already very much comparable to 
the ones at the low-index Pd(lOO) surface. Interestingly, 
the stability at the threefold (111) step sites is very simi- 
lar to the stability at the fourfold (100) terrace sites, and 
their energetic order depends on the local oxygen concen- 
tration. We speculate that this may lead to interesting 
ordering phenomena, on which we focus in ongoing work. 
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